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Abstract 

In this paper we investigate some aspects of the dynamics of a quantum system constituted by a 
two-level atom (atom 1) in interaction with a single mode cavity field. The field is considered to be 
weakly coupled to a third subsystem: a second two-level atom (atom 2), in a mixed state, which will 
play the role of an external environment. We analyze the influence of such a small environment on 
the dynamics of quantities such as the linear entropy (state purity) and the atomic dipole squeezing 
of atom 1, as well as the entanglement between atom 1 and the field. In order to study a general 
situation, we obtain the full analytical solution of the two atom Tavis-Cummings model in which 
the constituent atoms may be coupled with different strengths to the field and also have different 
frequency detunings. We show that nonclassical features associated to the main system may be 
significantly degraded by the action of the small environment, if atom 2 is resonantly coupled to the 
field. We also demonstrate that the nonclassical behaviour of the system may be restored if we detune 
the field from the transition frequency of atom 2, the environment. 


1 Introduction 

The evolution of quantum systems in contact with external environments has been the subject of inves¬ 
tigation for quite some time pQj . The environment is usually modeled as a (quantum) system having a 
large number of degrees of freedom, Viz., a reservoir, and number of methods, mainly perturbative, have 
been developed in order to study the influence of the surroundings on the behaviour of quantum systems 
[Tj. It is in general possible to derive an equation describing the evolution of the reduced density operator 
of the system of interest - the master equation, obtained by tracing over the degrees of freedom of the 
environment. Generally speaking, the coupling to an external environment has a detrimental action on 
the quantum properties of the system of interest, causing effects such as the loss of quantum coherence 
(decoherence) |2;[3J[3|. The models of environment mostly rely on the assumption of the existence of an 
ideal (large) reservoir, which naturally leads to irreversibility and relaxation features. Consider a single 
mode (cavity) field of the quantized field as being the system of interest. Cavity losses (dissipation) may 
be modeled via a reservoir constituted by a large collection of either independent field modes [I] or a 
beam of two-level atoms (5]. Both procedures lead basically to the same master equation for reduced 
density operator of the cavity field. Nevertheless, environments with other types of structure, not neces¬ 
sarily having a large number of degrees of freedom, may also occur in nature. Thus, one may think in a 
situation in which a system of interest is coupled to an uncontrollable surrounding environment having 
a small number of degrees of freedom. We could cite the interesting experiment (double photoionization 
of H 2 molecules) presented in reference | 6 ], where it is shown that even a single electron, constituting a 
“minimal environment” may strongly affect the interference fringes relative to a second electron (system 
of interest). The effects of small environments constituted by either a single harmonic oscillator or a single 
spin- 1/2 particle have also been addressed in [71 [5] and [S]- As a matter of fact, one of the smallest pos¬ 
sible environments could be constituted by a single two-level system, opposed to the model of a reservoir 
containing a large number of atoms [5j. In m it is studied a system made of a two-level system (atom 
1 ) coupled to an oscillator (field mode) which is itself in interaction with a minimal environment consti¬ 
tuted by a second two-level system (atom 2). This is formally the well-known two-atonr Tavis-Cummings 
model m but where an asymmetric partition of the system has been considered; the system of interest 
being constituted by atom 1 + field, while a partial trace is performed over the environment (atom 2 ). 
In [HE!, the discussion is restricted to the exact resonance case, i.e, the atom l( 2 )-field detunings being 
equal (Ai = A 2 = 0). Nevertheless the atoms are assumed to be coupled with different strengths to 
the field (Ai 7 ^ A 2 ); while atom 1 and the field are assumed to be strongly coupled, the field is weakly 
coupled to its environment. If the field is completely isolated from atom 2, i.e., A 2 = 0.0, we have the 
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basic model of quantum optical resonance known as the Jaynes-Cummings model [12] . a particular case 
of the Tavis-Cummings model for N = 1. The Jaynes-Cummings model has well known features; for 
instance, coherent Rabi oscillations of the atomic inversion as well as the linear entropy of the atom, 
S = l-Tr(p 2 a ).^\ if the field is initially prepared in a Fock state. Such a characteristic behavior may be 
useful to evaluate the interference of external systems on the regular evolution of a quantum system, and 
thus we could focus our attention on the reduced density operator of atom 1 , p a i, obtained from a partial 
trace over both the variables of atom 2 and the field. One of the conclusions of m is that the evolution 
of the linear entropy strongly depends on the degree of mixedness of atom 2 , meaning that even a small 
(although noisy) environment may significantly disturb the evolution of two coupled quantum systems 
(atom 1 + field). Here we would like to further investigate in which way the quantum features of the 
system of interest are affected due to its coupling to such a small environment. We analyze the evolution 
of non-classical effects like the atomic dipole squeezing Ha nn, as well the entanglement between atom 1 
and the field. In general non-classical effects are susceptible to external influences, and we would like to 
find to what extent an unwanted small environment with few degrees of freedom may affect the squeezing 
and entanglement properties of our system. We also address a more general situation, in which the detun¬ 
ings Ai, A 2 as well as the coupling strengths Ai,A 2 are different. This provides us additional flexibility. 
For instance, we would like to minimize the destructive effect of the environment, and one possible way 
of doing it is by controlling the detuning between the field and the environment itself. We note that 
the analytical solutions of the two-atonr Tavis-Cummings model found in the literature are restricted 
to particular cases. Namely, the atom-field coupling constants may be assumed to be equal (identical 
atoms) 15], or different (non-identical atoms), but still having the frequency of the field equal to the 
atomic transition frequencies mm- a step further is given in EH!, where it is presented a solution 
of the model for non-identical atoms but being both of them either on resonance or out of resonance 
with the field. Here we work out an exact analytical solution of the two-atom Tavis-Cummings model 
for distinct coupling constants, Ai,A 2 and arbitrary detunings, Ai, A 2 . We remark that the solution in 
this general case has a considerably more complicated structure than the solution for Ai = A 2 (equal 
detunings). 

Our paper is organized as follows. In section 2 we obtain the general analytical solution of the 
two-atom Tavis-Cummings model with arbitrary detunings and couplings. In section 3 we discuss the 
time evolution of the linear entropy, the atomic dipole squeezing of atom 1 as well as the atom 1 -field 
entanglement, considering atom 2 as a disturbance (environment). We also show how it is possible to 
restore quantum coherence, the atomic dipole squeezing and entanglement by controlling the frequency 
detuning between the field and the environment (atom 2). In section 4 we summarize our conclusions. 


2 Tavis-Cummings model: an analytical solution for different 
coupling constants and different detunings 

The two-atom Tavis-Cummings model is described by the following Hamiltonian (under the rotating wave 
approximation and making Ti = 1 ) 


H = —of + ~ a 2 + ua)a + Ai(acr j f + a)a 1 ) + X 2 (aa 2 + aV 2 ) 


(1) 


where a (a^) are the annihilation (creation) operators associated to the field mode, with frequency w; 
erf , cF 1- e < 7 - are the de Pauli operators relaive to the “i-th” atom, each atom having transition frequency 
uii . We may rewrite the Hamiltonian above in terms of the detunings Ai = u>\ — to and A 2 = 0 J 2 — uj, as 


+ a)a ) + ———< 7 ^ + —-erf + Xi(aa^ + a)a 1 ) + A 2 (aerf" + a) cr 2 ) . 


H = uj 

The Hamiltonian ([ 2 ]) may be split in two parts 

H = H 0 

where 

H 0 = u 


■H u 


cr$ +cr 2 t 
— 1 ——- + a'a 


( 2 ) 

( 3 ) 

( 4 ) 


2 Being p a the reduced atomic density operator. 
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represents a conserved quantity, and 




~^~ a i + _ 2^ <t 2 + Ai(acr+ + o’ 1 ’o'! ) + A 2 (oct^" + a/a 2 ) 


( 5 ) 


is the interaction part. 

The Schrodinger equation in the interaction representation reads 

i J t Hi (*)) = Hi Hi W) > ( 6 ) 

with Hj = exp (z H 0 t) H 1 exp (—z H 0 t), and Hj = H\. 

The Hamiltonian Hi induces transitions between the states |ei, e 2 , zi), |ei, g 2 , n + 1), \gi, e 2 , n + 1), 
\g±, g 2l n + 2), and therefore we may write the following ansatz 

Hi (*)) = C'l.nW |ei,e 2 ,n)+C' 2; „(i)|ei,52,« + 1) 

+ C 3 ,n (t) \gi, e 2 , zz + 1) + C 4 ; „ (t) |e/i, g 2 , n. + 2} . (7) 

After substituting the proposed solution ([7]) above in the Schrodinger equation, we obtain the correspond¬ 

ing set of coupled differential equations for the amplitudes Ci, n {t) 

iCl,n = ^- 2 -^ Cl ,n + A 2 \/zi + 1 C 2 , n + XlVn + 1 C3,n , 

iC 2y n = A 2 \/z? l + 1 Ci.n + ^^ C^n + AlVVl + 2 C4 jra , 

iCz,n = AlV^ + IC'l, n + (---^ C3,n + A 2 v'n + 2 C 4 jT j , 

iCi^n = X\\/n + 2<7 2 , n + A 2 Vzi + 2^C 4 iT i. (8) 

Because the exact analytical solutions of the equations ^ are rather lengthy, we have placed 
them in the appendix. 


3 Numerical results: quantum state purity, dipole squeezing 
and entanglement 

Now we would like to discuss the reduced dynamics of the quantum system, basically focusing on the 
properties associated to atom 1. In order to do so, we should first trace, the total density operator 
p(t), over the variables of atom 2 (the “environment”), obtaining the joint atom 1-field density operator, 
or p(t) a ij = Tr a 2 p{t). For a completely isolated atom (Jaynes Cummings model, [12]), we know that 
the linear entropy has a completely reversible behaviour, provided the field is initiallly prepared in a 
Fock state. Moreover, non-classical features such as squeezing m may also arise during the atom-field 
interaction. If we want to focus solely on the atomic properties, we should perform a further partial 
trace over the field variables, i.e., calculate the reduced density operator relative to atom 1, p(t ) a i = 
Trf [p{t)aij]- Here we are going to assume initial conditions of the form: p (0) = p a i (0 )®pf (0 )®p a 2 (0), 
with 

Pal (0) = Hi) Hi\ , pf (0) = \N) {N\ , p a2 (0) =p|e 2 ) (e 2 | + (1 - p) \g 2 ) (g 2 \ ■ (9) 

In other words, having atom 1 prepared in the state = cos (§) |< 7 i) + sin (|) e 1 ^ |ei) and the field 
prepared in a Fock state | N), with N being either N = 1 or N = 0. Given that atom 2 is treated here as 
a noisy “environment” perturbing the atom 1-field system, we consider it initially in a statistical mixture 
of its ground and excited states, p a2 (0); here we are going to address the case of a maximally mixed 
environment, or p = 1/2. We may now discuss some of the quantum dynamical features of the system. 
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3.1 Reduced dynamics of atom 1: linear entropy evolution 


The linear entropy is used to quantify the degree of purity of a quantum state. In the case of atom 1, it 
may be written as 

S{t) = l- Tr al (p 2 al ) = 2 [a - a 2 - | 7 | 2 j , (10) 


where 


and 


a (t) = sin 2 


+ cos I - ) p 


l '3,1 °4,1 


c: 


(hi) 


3,0 


+ 


c: 


(in) 


4,0 


+ (1 ~P) 


+ (1 ~P) 


a 


3,0 


a 


(in 

4,0 


c. 


(IV) 

3,-1 


r (IV) 

°4,-l 


7 W = 


„-i^, sin (9) 


r (iy r (in) r (iy r (ui) 

°1,1 °3,0 °2,1 °4,0 


are the populations and coherences (respectively) of the reduced density operator of atom 1 


( 11 ) 


( 12 ) 


Pa\ ( t ) = a\g 1 ) (pi| + (1 - a) |ei) <e x | + 7 I 51 ) <e x | + 7 * |ei) {gi\ . (13) 

Here the superscript i, (i = I, II, III, IV) in the amplitudes C^ n in equations (JTTJ) and |l 2 | above are 
related to the initial conditions. The amplitudes Cj iTl (t) may be written as a linear combination of the 
initial conditions as Cy n (t) = m _ 1 Aj m Cj (0). The notation employed is such that C^J, for instance, 

indicates that C\ (0) = 1 and the remaining coefficients are zero; the superscript (II) in c\\ ^ indicates 
that C 2 ( 0 ) = 1 , and the remaining coefficients are zero and so on. 

We would like now to analyze the effect of A 2 (atom 2-field detuning) on the time evolution of 
the linear entropy of atom 1. In figure |l]) we have plots of the linear entropy of atom 1 as a function of 
time considering different values of A 2 . The field is initially prepared in the one photon Fock state |1), 
atom 1 in its excited state |ei), and atom 2 in a maximally mixed state, p a 2(0) = | (\g 2 ) (<72! + 162) (e 2 1 ) ■ 

As seen in figure the linear entropy is a periodic function of time (Rabi oscillations) if 

atom 2 is decoupled from the held, i.e., for A 2 = 0.0. If the interaction between atom 2 and the held is 
turned on (f A 2 = 0 . 1 ), the evolution of the linear entropy becomes very irregular, as shown in hgure @>)> 
which characterizes a destructive effect due to an unwanted coupling to a noisy sub-system. Nevertheless, 
it is possible to restore periodicity by controlling the atom 2-field frequency detuning. If A 2 is increased, 
we expect less influence of atom 2 over the dynamics of the system of interest. In fact, as shown in hgure 
0 ) and hgure ©). periodicity may be re-established for a sufficiently large A 2 . We remark that atom 
1 is kept on resonance with the held (Ai = 0.0) as well as strongly coupled to it (Ai = 1.0). 


3.2 Reduced dynamics of atom 1: atomic dipole squeezing 

Atomic dipole squeezing is a non-classical effect that may arise in the Jaynes-Cummings model [13] . 
An atom is said to be dipole squeezed if the quantum fluctuations of the atomic dipole are below the 
fundamental limit imposed by the Heisenberg inequalities. The components of the (slowly varying) atomic 
dipole operator (for atom 1) may be written as [13] 


a x = ate iUlt + aie~ iuit , 

o y = ^ (a+e iuit - . (14) 

The operators above do not commute ([<r x ,a y ] = 2ia z ), and thus they should obey the Heisenberg 
inequality (A cr x ) (A a y ) > |(cr z )|. Atomic dipole squeezing is verihed if 

(Acr x , y ) 2 < I ( 07 ) I , (15) 
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where (A cr x ) 2 = 1 — 4(i?e(cr 1 )e IWl *) 2 and (Acr y ) 2 = 1 — 4(/to(ct 1 )e IWlt ) 2 . The conditions for 
squeezing may be written in terms of the functions si e S 2 (indexes of squeezing) as 


si = 


1-4 (Re(ai)e- iu>lt y 

KOI 


<1 or s 2 = 


1 - 4:(lm(ai)e- iu>lt y 

10)1 


< 1, 


(16) 


where oji is the frequency of atom 1 , and (cr z ) is its corresponding atomic inversion. Here (<r Q ) = 
Tn ( (T a p a i )• 

We may rewrite the indexes of squeezing, S! and s 2 , in terms of the functions a(t) and 7 (t) 

above, or 


si 


s 2 


1 — 4 (Re[ 7 ] cos (coit) — Im[ 7 ] sin (uqf )) 2 
11 2ck (£) | 

1 — 4 (i?e[ 7 ] sin (uq t) + Im[ 7 ] cos (uqf )) 2 

|l-2a(i)| 


(17) 


Now we choose specific initial conditions which allow dipole squeezing in the case of having atom 
2 completely decoupled from the field (A 2 = 0.0), which corresponds to the Jaynes-Cummings model. 
Atom 1 is assumed to be prepared in the superposition state |0i(0)) = cos (0.6) \g\) + sin (0.6) |ei); atom 
2 in the maximally mixed state, p a 2 (0) = \ (|<? 2 ) (g 2 | + |e 2 ) (e 2 |) and the field in the one photon Fock 
state, |1). In figure ([ 2 ]) we have plotted the dipole squeezing index Si relative to atom 1 as a function 
of time. We note that in the absence of the “environment” (atom 2), dipole squeezing occurs for a few 
narrow intervals of time [figure (J2^,)]. However, if atom 2 is weakly coupled to the field, e.g., A 2 = 0.1, 
the dipole squeezing is inhibited due to the action of the environment [see figure ©3)]- Nevertheless, 
similarly to what we have seen in the previous subsection, dipole squeezing may be restored for a large 
enough detuning between atom 2 and the field, as shown in figures ©;) and #)■ 


3.3 Atom 1-field entanglement 

We may quantify the entanglement between atom 1 and the field by evaluating the negativity m , an 
entanglement measure especially convenient in this case. Now we assume the following initial conditions: 
the field initially prepared in the vacuum state 10 ), atom 1 in its excited state |ei), and atom 2 in a 
maximally mixed state, p a 2 (0) = ^ (|p 2 ) (<? 2 | + |e 2 ) (e 2 |). The negativity may be calculated from the 
time-dependent atom 1 -field reduced density operator 


Pal-/ (t) = P11 |ffi, 0) (gi, 0| + P 22 |si, 1) (<71, 1| + P33 Iflij 2) (<71, 2| + 

+P44 |eii 0) (ei, 0| + P55 |ei, 1) (ei, 1| + p 2 4 |pi, 1) (ei, 0| + 
+P24 l e i> 0 ) {dh 1 | + P35 Iflii 2 ) (ei, 1 | + P3.5 |ei, 1 ) (<71, 2| . 


The relevant matrix elements are 


1 


P11 = 


a 


(a) 


3,-1 


P 22 = 


a 


+ a 


Mi) 

y 4,-l 


P33 = 


fi(I) 

l '4,0 


P44 = 


u l,0 


c: 


(a) 


P55 = 


c. 


p 35 = 


K„>y 


c 


(I) 

4,0 J 


k [KK 


c. 


(i) 


P 24 — 2 LV Wl .°) w 3,0 
Thus the negativity M may be expressed as 


(cH)'cS 


-V = 2 E 


i= 1 
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°1 — 2 ( 'p 22 ~ \!P 22 + 4 |P 35 1 2 ^ , — 7 , ^Pll + P 55 - \/(pil - P55 ) 2 + 4 |P 24 1 2 

We have inserted a factor of 2 so that the maximum of entanglement corresponds to 1. We 
note that although the negativity is somehow related to the linear entropy, they may not have the same 
behaviour, given that in this case the reduced atom 1 -field density operator is obtained from a mixed 
state. In figure ([3]) we have plotted the negativity J\f as a function of time. For the initial conditions 
chosen here, the system has a simple evolution, i.e., a oscillatory behaviour as shown in figure § 0 , 
provided the system is isolated (decoupled environment). However, if atom 2 (environment) is resonantly 
coupled to the field, the entanglement is severely damped [see figure (§>)], similarly to what happens to 
the linear entropy and atomic dipole squeezing. Again, entanglement may be restored if A 2 is increased, 
as shown in figures (|3jc) and ([3jl). 

4 Conclusions 

We have presented a study of the dynamics of a bipartite quantum system (atom 1 -field) in which the 
field is weakly coupled (A 2 « Ai/10) to a small environment (atom 2). A single atom corresponds to the 
smallest possible environment considering the atomic beam model for a reservoir [5]. In the many atoms 
model a Markovian master equation is obtained via a perturbative method, which leads to irreversible 
field damping and decoherence mm- However, in the “single atom environment” model an exact, non- 
perturbative solution is allowed, and as found in [ 10 ] (restricted to the resonant case), the linear entropy 
of the quantum system undergoes an irreversible-like evolution which is very sensitive to the mixedness of 
the small environment. Here we have considered a more general situation, by analytically solving the two- 
atom Tavis-Cummings model for non-identical atoms, with different coupling constants, (Ai 7 ^ A 2 ) and 
different detunings (Ai 7 ^ A 2 ). Besides, we have also investigated the evolution of important quantities of 
the system related to non-classical behaviour, such as the dipole squeezing of atom 1 and the negativity, 
which quantifies the atom 1-field quantum entanglement. The small environment (atom 2), has been 
assumed to be initially in a maximally mixed state. Naturally in this simple model, contrarily to what 
happens in the case of a many atom reservoir, there is no irreversible loss of coherence. Instead, due to the 
incommensurate frequencies characteristic of the model, the resulting dynamics is quasi-periodic, and we 
should expect quasi-recurrences of the considered physical quantities at longer time-scales. Nevertheless, 
due to the incoming “thermal noise” from the environment, the atomic dipole squeezing as well as the 
atom 1 -field entanglement are considerably degraded (in the resonant case, A 2 = 0 . 0 ), specially if one 
is restricted to shorter time scales (up to a few Rabi cycles, as shown in the figures). In particular, we 
have verified that dipole squeezing is completely suppressed. Squeezing only exists if the index defined 
in equation ( [l7| becomes less than one [see figure and it occurs at relatively narrow time intervals. 

We have also obtained a full analytical solution of the two-atom Tavis-Cummings model for 
non-identical atoms, which allowed us to assess the effect of the atom 2 -field detuning on the dynamics of 
the system. Indeed there is a competition between the environment-field coupling and the environment- 
field detuning, but although we were able to determine explicitly the dependence of the quantum state of 
the system on the parameters (Ai, A 2 , A lt A 2 ), the lengthy expressions hindered a more detailed analysis. 
Yet, it is possible to have some degree of control over the dynamics of the system and circumvent the 
destructive effect of the environment. By increasing the atom 2-field detuning A 2 , there will be an effective 
decoupling of the small environment, and we expect a restoration of the non-classical properties such as 
squeezing and entanglement. For instance, we have found that it is sufficient to have a (environment-field) 
detuning around A 2 = 1.0 in order to recover the atomic dipole squeezing property of atom 1. However, 
larger values of detuning would be required to restore the atom-field entanglement, as it is clearly seen 
in figure §]• We shall also point out that, as discussed in m, any additional disturbances combined to 
the action of the small environment will turn out to be very detrimental for quantum coherence. 
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Appendix. Solution of the two-atom Tavis-Cummings model 


Here we present the full analytical solution of the two-atom Tavis-Cummings model. In order to solve 
the system of coupled differential equations ([s]) , we have employed the Laplace transform method, and 
the set of differential equations is transformed to the following set of algebraic equations 


^---b z Ci — \ 2 y/ n + 1 Ci — Ai y/n + 1 C 3 

—Ai \/n+ 1 Ci + + is^j C 2 - Ai \Jn + 2 C 4 

-Ai A + 1 C-\ + 2 + is^j C 3 - \ 2 Vn + 2C 4 

—Ai y/n + 2 C 2 — A2 Vn -I- 2 C 3 + (-—-1- i s j C 4 


ici (0) 
ic 2 (0) 
ic 3 (0) 
i C4 (0). 


For simplicity we have denoted Cj (s) = Cj in the equations above. The subsequent steps 
involve the solution of polynomials up to fourth degree, and after some involved calculations we obtain 
the full solution. 

The diagonal terms are 


e -5 ((n/X-vTT) ( 18 (A + z) DiJ ± (-2zfci + 3 vA - 3Ai) - 4320/ - Ai.J'f + (Ay - 3i/X + 3 Ad) 3 ) 

AT (-4 AT AT + 4X + y - y) 

e-HCx+VYl) ^ 18 ^3 + iDi y U ki + 3i(yx + at) ) + 432£>3 + i ( 4 Jf + (2ki + 3z {^X + AT) ) 3 ) ) 

AT (4vAAT + 4X + y - y) 

glA-A (i 8 (A + i) Di Ji (2ifci + 3 VT - 3 At) + 4320? + 4zJ? + (-2 iki - 3 AT + 3 A 2 ) 3 ) 

AT (-4VAAT + 4 x-y 1 + y) 


e 4<(vx-h/A | 

()l8 (1 - zA) Or Ji | 

(Ai — 3z 

(AT + At) 

) - 4320? - z 

(4 J? +1 

^2fci — 3z | 

(AT + AT) 

) 3 ))) 



AT ( 4 ATAT + 4X - y + y) 



An — 


54 


e -|t(v^-v/Yi) ( 18 (A + z) 0 2 J 2 (—2zfc 2 + 3y/X — 3AT) - 4320? - 4i,/| + (Afc 2 - 3vA + 3 AT) 3 ) 

AT (-4vAAT + 4X + y - y 2 ) 
e -H'/y+Cn) (i 8 (A + z) iD 2 J 2 (2k 2 + 3 i [yx + AT)) + 4320? + z (4 Jf + (Tfc 2 + 3z (VX + AT)) 3 )) 

AT (4vAAT + 4X + y - y ) 


e 5t(v^-v^) | 

(l8(A + i) AA 1 

(2zfc 2 + 3\/X — 3 AT) 

| + 4320? + 4z J? + | 

(-2zfc 2 - 3y/X + 3 AT) 

3 ) 

AT (-4vAAT + 4X - y + y) 

1 


e i t (Vx+VK) | 

(is (1 - zA) o 2 j 2 1 

(2fc 2 — 3z 

(at+at) 

) - 4320? - i 

(44 +1 

(2fc 2 — 3z | 

[y/X + AT) 

) 3 ))1 



AT ( 4 Vx AT + 4X - y + y) 



A22 "54 
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x r (l 8 (v^ + i) D 3 J 3 (- 2 ik 3 + 3VX - 3^) - 432L>f - 4 iJ§ + ( 2 ik 3 - 3^ + 3-v/Yi) 3 ) 

54 \ VV (-4 VxVV + AX + Y, - Y 2 ) 

e - i t (Vx+VYI) ^ lg (^3 + ^ i£>3 J 3 (2fc 3 + 3* (Vx + y/YiJ ) + 432D| + i ( 4 J| + (2 k 3 + 3* (v 7 ^ + y/Yi) ) 3 ) ) 

VYi (4\/XVYi + 4 X + Y-l- Y 2 ) 

e H^X-VYi) ( 18 (yg + j) D 3 J 3 ( 2 ik 3 + 3VX - 3v 7 ^) + 432D| + 4i Jf + (-2ifc 3 - 3\/X + 3\/Y^ 3 ) 

V%. (-4^^ + 4X - Yi + F 2 ) 

e |t(v / x+v / ^) ( 18 _ y 3 ) D 3 J 3 (2k 3 - 3 i (VX + v 7 ^) ) - 432Df - i (4Jf + (2fc 3 - 3i (v 7 ^ + v 7 ^) ) 3 ) ) j 

VYi (4 VXVY 2 +AX-Y 1 + Y 2 ) J 


1 r e -5*('/5f-v^r) ( 18 + j) £) 4 J 4 (_2jfc 4 + 3\/X - 3^) - 432D| - 4* Jf + (2ifc 4 - 3^ + 3^) 3 ) 

54 \ VYi (-4-VxVYi + 4X + Y- y 2 ) 

e -H n/x+v^T) ( 18 y/3 + ^ iDi J 4 ^2fc 4 + 3i (v/X + v 7 ^)) + 432DI + i (d,/f + (2fc 4 + 3* (VX + VYi )) 3 )) 

v 7 ^ ( 4 ^^ + 4x + y - y 2 ) 

e | t(Vx-VY2) ( 18 ^y/3 + ;) d 4 j 4 ^2ifc 4 + 3VX - 3v 7 ^) + 432£>| + 4*Jf + (-2ifc 4 - 3^ + 3-v/yj) 3 ) 

VYi (-4 VxVYz + 4X-Y 1 + y 2 ) 

e |t(v / x+v^) ^ 18 ^ _ ix /3 ) d 4 J 4 ^2fc 4 - 3* (v 7 ^ + V 7 ^)) - 432D| - i (4Jf + (' 2k± - 3 i (VX + VY^j ) 3 )) 1 

V% (4VXVY2 + 4X-Y 1 + y 2 ) j 

and the non-diagonal terms read 


14 — — 4 AiA 2 y 7 (rT-l-T) 0 n-P 2 ) 

+ 


j (Vn-Vx) e-V(^-Yn) (yx + VYi) e~V(YX+VYi) 

\ v 7 ! 7 ! [~4VxVYi + 4X + Y- y 2 ) v 7 ^ (dv^v 7 ^ + 4X + y - y 2 ) 

(VX + VY^eHYx^) (^VYi - Vx^eViYx-VY) j 

V 7 ^ (dv^v 7 ^ + 4X- Yi + y 2 ) v 7 ^ (-4 VxVYi. + 4X - Yi + Y 2 ) J ’ 


^2 3 — — 


2e it(^+v^) ^ ( Al + a 2 ) AiA 2 + (2 n + 3)y/x + (2 n + 3) v 7 y 


2 + 4X - Y 4 + Y 2 ) 
e .''(v'A'VV:) ^ 2 - ( Al + A2 ) A : A 2 - 2(2n + 3)VX - 2(2n + 3)-^) 


1 + 4X + Y 4 - Y 2 



git(v'x-v'Yh) ^2 i (Ai + A 2 ) AiA 2 + (4 n + 6 )VX - 2(2 n + 3 )y/Y 2 s j 
v 7 ^ (-4 yfXy/%. + 4X — Yi+ y 2 ) 
g-i^v^-v^T) (j 4n + 6)a /x _ 2 ((2 n + 3 )>/?[ + i (A x + A 2 ) A^)) 
y/Y x (-4 ^Xy/Y x +4 X + Yi- y 2 ) 


A 13 = —iXiVn+ 1 


+ 


r (g+(VX-VY 2 ) (- 2 iA 2 + VX-VY 2 ^ 

\ VY> (-Ay/Xy/Yi + 4X — Y 1 + Y 2 ) 

git(yx+v^) ( G + + ^/yfj 2iA 2 + v 7 ^ + v 7 ^)) 

v 7 ^ (4^^ + 4X - n + y 2 ) 

g-|t(vX-v^T) (g+ (v^-v 7 ^) (2iA 2 + >/X-VTT)) 

v 7 ^ 4 v 7 X v 7 4T + 4X + y - y 2 ) 
e - 3 t ('/A+v / YT) ^_ G _ ((2*A 2 + yfx + Vy) ) | 
V 7 ^ ( 4 ^^ + 4X + y - y 2 ) 


^4i2 = — iy/ n + 1A 2 


+ 


gi t(Vx-SK) (q _ (yx - y/Y^j (-2iA, + v 7 ^ - ) 

v 7 ^ (-4VX^ + 4X-Y 1 + y 2 ) 
gHvX+v^) G + (yfX + v 7 ^) (-2iAi + VX + v 7 ^ 
v 7 ^ (4^X^ + 4X-Yi + y 2 ) 


g“ ^(v 7 ^ - (- G + (%/X - v 7 ^) (2iAi + \/X - v 7 ^) ) 
^/Y x 4v / X v / y( + 4X + Yi - y 2 ) 
e -*t(VxW57) ( G _ ^/x + ^ (2iA 1+ VX + ^)) | 
V 7 ^ ( 4 ^^ + 4X + y - y 2 ) 


A 2 4 = —iy/n + 2Ai 


[ r. •M/- v Vv,) (# + (\/X-v 7 ^) (-2 zA 2 + v 7 X- v 7 ^)) 

[ VYi (~4y/Xy/Yi + 4X + Yi - y 2 ) 

e -i t (vX+yiT) _ (Vx +v 7 ^) (-2iA 2+1 /X + yy)) 
V 7 ^ ( 4 ^^ + 4X + Y 1 - y 2 ) 
e HYx+VY) (y + (yx + ^Y^j ( 2 ; a 2 + \/x + v /^)) 
V 7 ^ ( 4 ^^ + 4X - Y, + y 2 ) 
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e i t (vW-vW) ( 

ff +1 

(Vx-VY) 

(2iA 2 + v / X-v / ^)) 

(-4\/X\/y + 4X - Yi + y 2 ) 

1 


^•34 — — Wti H - 2A 2 


j e-hKVx-YY) | 

(-ff + 1 

(Vx-VY) 

(- 2 iA 1 + Vx-VY)) 

| v /Y (-4VXVY +4 X + Y 1 - F 2 ) 

1 


^ >/X + >/TT (—2zA 


+ 4X + y - y 2 ) 


c ',t(YX s/Y,) ^ 

H — | 

(Vx- v 7 ^) 

(2iAr + y/X- VY)) 

JY 2 (-AjXyfYi + 4X-Y 1 + y 2 ) 

1 


,i*(^WW) (_ H + fyx + VY) ( 2 *Ai + \/x + v^)) 


l 2 ^‘±V XV V -r 2 + 4X — Y 1 +Y 2 
In the expressions above, the coefficients X and Y\ t 2 are given by 

2 A 2 1 / 3 (A 2 + 12 F) r 


X = — 


Y 12 = -2 A-X± 


3/3 3x2V3 

2iB 


Vx' 


where 


r = 


2xl 3 - 27B 2 - 72AF + \J -4 (A 2 + 12F) 3 + (2 A 3 - 27B 2 - 72AF) 2 


1/3 


and 


A — - (A 2 + A 2 + 2(2n + 3) (Aj + A 2 )) 


b = (a 1 x 2 2 + a 2 \ 2 1 ), 


F = w^ A ^ A2 > 


(A 2 - A^) 2 + 4 (3 + 2 n) (Af - A 2 ) + 16 ( n 2 + 3?i + 2) (A 2 - A^) 2 } 


We have also that 

G = A 2 — A 2 + 4 (n + 2) (Aj — A 2 ) , 
H = A 2 — A 2 + 4 (n + 1) (A 2 — Aj), 


A = 


(-1) 5/6 (3a ;: + fc 2 ) 
3 Ji 


i = 1, 2, 3, 4, 
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and 


Ji = (-27 bi + 9 (nki + 2 kf - 3 iV3y/4a% - 27 b 2 + 18 a^h + a 2 k 2 + 46 ,' , * = 1 , 2 , 3 , 4 , 


being 


ai = a4 = 


(Ai - A 2 ) 2 + 4 (n + 2) (A? + \ 2 2 ) 


a 2 = a 3 = 


(Ai + A 2 ) 2 + 4 ((n + 1) Af + (n + 2) A 2 ) 


h = -6 4 = - (Al g A2) [A? - A| + 4 (n + 2) (A? - A 2 )] 


b -2 = -bs = - (Al + Aa) [a; - A 2 + 4 ((n + 1) A 2 - (n + 2) A 2 )] 


fci = — AT4 = — 


(Ai + A 2 ) 


and 


k 2 = -fc 3 = - 


(Ai — A 2 ) 
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Figure 1: Linear entropy of atom 1 as a function of time. In (a) the field is not coupled to the environment, 
i.e., A2 = 0.0 (Jaynes Cummings model). For the field coupled to the environment, with A2 = 0.1 and 
(b) A 2 = 0.0; (c) A 2 = 1.0, and (d) A 2 = 5.0. In all cases Ai = 1.0, and the initial state of the system is 
a tripartite product state with |-f/> a i) = |ei), p a 2 = \ (I 52 ) ( 32 I + |e 2 ) (e 2 1 ) and | t/jf) = |1). 
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Figure 2: Squeezing index si as a function of time. In (a) the field is not coupled to the environment, 
i.e., A2 = 0.0 (Jaynes Cummings model). For the field coupled to the environment, with A2 = 0.1 and 
(b) A 2 = 0.0; (c) A 2 = 1.0, and (d) A 2 = 5.0. In all cases Ai = 1.0, and the initial state of the system 
is a tripartite product state with|?/> a i) = cos (0.6) \g\) + sin (0.6) |ei), p a 2 = * ( \g 2 ) (52I + \e 2 ) (e21) and 

l^/> = I 1 )- 
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Figure 3: Negativity J\f used to quantify the entanglement between atom 1 and the field as a function 
of time. In (a) the field is not coupled to the environment, i.e., A2 = 0.0 (Jaynes Cummings model). For 
the field coupled to the environment, with A2 = 0.2 and (b) A 2 = 0.0; (c) A 2 = 1.0, and (d) A 2 = 5.0. 
In all cases Ai = 1.0, and the initial state of the system is a tripartite product state with \ipai) = |ei), 
Pai = 5 (|S2> (521 + |e 2 ) (e 2 |) and \ip f ) = |0). 


14 






[7] G.-L. Ingold, P. Hanggi and P. Talkner, Phys. Rev. E 79 (2009) 061105 

[8] A. Vidiella-Barranco, Physica A 459, 78 (2016) 

[9] M. Castagnino, S. Fortin and O. Lombardi, Mod. Phys. Lett. A 17, 1431 (2010) 

[10] A. Vidiella-Barranco, Physica A 402, 209 (2014) 

[11] M. Tavis and F.W. Cummings, Phys. Rev. 170, 379 (1968) 

[12] E.T. Jaynes, F.W. Cummmings, Proc. IEEE 51, 89 (1963); B.W. Shore, P.L. Knight, J. Mod. 
Optics 40, 1195 (1993) 

[13] D.F. Walls, P. Zoller P, Phys. Rev. Lett. 47, 709 (1981); X.S. Li, D.L. Lin, T.F. George, Phys. 
Rev. A 40, 2504 (1989) 

[14] T. Maqbool, M.S.K. Razmi, J. Opt. Soc. Amer. B 10, 112 (1993) 

[15] Z. Deng, Opt. Comm. 54, 222 (1985) 

[16] S. Mahmood, M.S. Zubairy Phys. Rev. A 35, 425 (1987) 

[17] H.T. Dung, N.D. Huyen, J. Mod. Optics 41, 453 (1994) 

[18] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, Phys. Rev. A 58, 883 (1998) 


15 



